Astronomy & Astrophysics manuscript no. Ietter'hdl90073'v5'ee 


©ESO 2013 


January 10,2013 





Letter to the Editor 



The dramatic change of the fossil magnetic field of HD 190073: 
evidence of the birth of the convective core in a Herbig star ? * 

E. Alecian 1 , C. Neiner 1 , S. Mathis 2 ' 1 , C. Catala 1 , O. Kochukhov 3 , J. Landstreet 4 - 5 , and the MiMeS collaboration 



CO 

o 

CN 
C 

a 

ON 



6 



> 
o 

00 



O 

CO 



X 



LESIA-Observatoire de Paris, CNRS, UPMC Univ., Univ. Paris-Diderot, 5 place Jules Janssen, F-92195 Meudon Principal Cedex, 
France, e-mail: evelyne.alecian@obspm.fr 

Laboratoire AIM, CEA/DSM, CNRS, Universite Paris Diderot, IRFU/SAp Centre de Saclay, F-91191 Gif-sur-Yvette, France 
Department of Physics and Astronomy, Uppsala University, Box 516, SE-751 20 Uppsala, Sweden 
Armagh Observatory, College Hill, Armagh, BT61 9DG, Northern Ireland, UK 

Department of Physics and Astronomy, The University of Western Ontario, London, Ontario, N6A 3K7, Canada 



Received September 15, 1996; accepted March 16, 1997 



ABSTRACT 



In the context of the ESPaDOnS and Narval spectropolarimetric surveys of Herbig Ae/Be stars, we discovered and then monitored the 
magnetic field of HD 190073 over more than four years, from 2004 to 2009. Our observations all displayed similar Zeeman signatures 
in the Stokes V spectra, indicating that HD 190073 hosted an aligned dipole, stable over many years, consistent with a fossil origin. 
We obtained new observations of the star in 2011 and 2012 and detected clear variations of the Zeeman signature on timescales of 
days to weeks, indicating that the configuration of its field has changed between 2009 and 2011. Such a sudden change of external 
structure of a fossil field has never previously been observed in any intermediate or high-mass star. HD 190073 is an almost entirely 
radiative pre-main sequence star, probably hosting a growing convective core. We propose that this dramatic change is the result of 
the interaction between the fossil field and the ignition of a dynamo field generated in the newly-born convective core. 
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1. Introduction 

A few percent of main-sequence and pre-main sequence (MS) 
stars of intermediate mass host strong magnetic fields, of sim- 
ple (topologicall y dipolar) configuration, sta ble over many years 
or decades (e.g. Donati & Landstreet 2009). Recent studies of 
Herb ig Ae/Be stars (lAlecian et al.l I2012I) as well as of giants 
dAuriere et al.ll2.008l I2009I) have demonstrated a fossil link be- 
tween the magnetic stars at different evolutionary stages: from 
the pre-MS to the post-MS. The physical processes at the origin 
of these fields are starting to be understood. Due to the presence 
of these fields in very young stars, and to the lack of charac- 
teristics that could be linked to a dynamo, these fields are be- 
lieved either to be remnants of the Galactic fields that would 
have been swept up by the contracting core, or to have been gen- 
erated by a dynamo when part of the star was convective dur- 
ing the very early phases of it s life. This is the fossil theory. 
Brait hwaite & Nord lund (2006) have performed numerical sim- 
ulations showing that an arbitrary initial field quickly evolves 
into a relaxed stable configuration (see also iDuez & M athis 
l2010t IDuez et alJ l20 10) and, once formed, evolve on timescales 
longer than the MS lifetime of a star. The resulting config- 
urations are similar to those observed on the surface of the 
intermediate-mass stars. One of the most poorly understood as- 



* Based on observations collected at the Canada-France-Hawaii 
Telescope (CFHT) which is operated by the National Research Council 
of Canada, the Institut National des Sciences de l'Univers (INSU) of 
the Centre National de la Recherche Scientifique (CNRS) of France and 
the University of Hawaii, at the Observatoire du Pic du Midi (France), 
operated by INSU, and at the European Southern Observatory, Chile 
(Program ID 187.D-09 17) 



pects of the fossil scenario is the interplay between the dynamo 
magnetic fields arising in the convective envelope at the begin- 
ning of the pre-MS phase, and in the convective core at the end 
of the pre-MS, and the fossil magnetic field residing inside the 
radiative layers of the star. 

HD 190073 is an Herbig Ae star with an effective temper- 
ature of 9250 K dAcke & Waelkensl |2004|) . Its magnetic field 
was discovered in 2 005 as part of a sur vey of magnetic fields in 
Herbig Ae/Be stars (ICatala et al.ll2007l) . We first observed it for 
4 consecutive days, then, once every few months for more than 
four years, with ESPaDOnS or Narval, to determine the stability 
of the field over time. In contrast to most magnetic V signatures 
observed in intermediate and high-mass stars, the Zeeman signa- 
ture of HD 190073 was always similar from one observation to 
another. No variability was detected. In other stars the observed 
variations are explained by the oblique rotator model, desc ribing 
a rota ti ng star hosti n g an i nclined dipolar magnetic field (Stibbs 
119501) . ICatala et al.l (120071) proposed three scenarios to explain 
the absence of rotational modulation in HD 190073: (i) the star 
is seen pole-on, (ii) the magnetic and rotation axes are aligned, 
or (iii) the rotation period of the star is very long. 

The stability of magnetic fields in Herbig Ae/Be stars, in- 
cluding HD 190073, is well establish ed, and is a very strong 
argum ent in favour of the fossil theory (Braithwaite & Nordlund 
2006). However, in 2011 and 2012 we obtained a new series of 
observations (about two years after the end of the first series), us- 
ing HARPSpol and Narval, and to our amazement we detected 
sttong variations on timescales of few days or weeks, clearly in- 
dicating an abrupt change of the magnetic configuration at the 
surface of HD 190073. This phenomenon has never been ob- 
served in any magnetic intermediate or high mass star. In this 
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Fig. 1. LSD V profiles (black) of HD 190073 ordered chronologically from the top left to bottom right. The numbers are the integer 
part of the HJD. Both 6176 profiles were obtained on the same day, 3.4 h apart. An error bar is superimposed on each pixel. The 
red profile overplotted on each observation is the average profile for the 2005-2009 data set. The letters refer to the instrument (E: 
ESPaDOnS, N: Narval, H: HARPSpol). 



letter we report the observations (Sect. 2), describe the analysis 
of the polarised spectra (Sect. 3), and discuss the possible origins 
of this unprecedented phenomenon (Sect. 4). 

2. Observations 

We have obtained two sets of observations. The first consists 
of 18 observations that were acquired between May 2005 and 
September 2009 with the high-resolution spectropolarimeters 
ESPaDOnS and Narval, installed at the Canada-France-Hawaii 
Telescope (CFHT, Hawaii), and at the Telescope Bernad Lyot 
(Pic du Midi, France), respectively. These observations were ac- 
quired in the course of PI programs dedicated to the study of 
the configuration and stability of the magnetic fields in Herbig 
Ae/Be sta rs. Some of t hese o bservations have already been pub- 
lished in Cata la et al.l (120071) . The second set of high-resolution 
spectropolarimetric observations consists of one HARPSpol 
(3.6m-ESO Telescope, Chile) spectrum obtained in May 2011, 
and 46 Narval spectra acquired between July 2011 and August 
2012. These observations were obtained in the framework of the 
MiMeS (Magnetism in Massive Stars) Large Programs, which 
are aimed at establishing the magn etic properties of massive 
stars at different evolutionary stages dWade et al.l l2009). 

All ESPaDOnS and Narval data were reduced as de- 

ol Libre 

2006). The HARPSpol observation was 
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scribed b' 
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red uced using a modi fi ed ve rsion of the REDU CE package 
of Piskunov & Valentil d2002l) . as described by iNeiner et al.l 
(120121) . In order to increase the signal-to-noise ratio (SNR) 



of our data we use d the Least-Squares Deconvolution method 
dPonati et al. 1997L LSD) wit h the complete Kurucz mask de- 
scribed in Catala et al .] (120071) . i.e. the mask containing almost 
all lines predicted by the model. This method provides us with 
the intensity / and circularly polarised Stokes V mean profiles of 
the spectra. These V profiles a re plotted in Fig. [T] Examples of 
LSD / profiles can be found in lCatala et all ( 120071) . The SNR in 
the LSD V profiles vary from 7000 to 22 000 in the two data sets. 
The longitu dinal magnetic field st rengths (Be ) were computed as 
described in lAlecian et al. I d2012h . and vary between -35 G and 
+55 G. 

The intensity spectra and LSD profiles of the second data 
set have similar c haracteristics than those of the first set, al- 
ready described bv lCatala et al.l(l2007l) . In both sets the spectrum 
is strongly contaminated with circumstellar emission, and this 
emission changes with time in a similar way. The photospheric 
(/) spectral lines of the two data sets have the same shapes, 
depths, and weak broadening. Unlike the intensity spectra, the 
circularly polarised spectra are quite different between the two 
data sets, as detailed and discussed below. 

3. The variability of the polarised spectra 

In Fig. Q] we plot (in black) the Stokes V LSD profiles obtained 
between 2005 and 2012. We have overplotted on each profile the 
average of all profiles obtained between 2005 and 2009 (in red), 
i.e. all profiles of the first observing set, during which no signifi- 
cant variations were observed. It is clear from Fig. Q~]that all pro- 
files of the first set are indeed similar to the average. However, in 
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Fig. 2. Average (top) and standard deviation (bottom) of the LSD 
V profiles of the first (black, 2005-2009) and second (red, 201 1- 
2012) data sets. 



201 1 and 2012, many profiles look different: sometimes the sig- 
nature is smaller in amplitude (e.g. HJD 2455753 and 2455790), 
sometimes there is no signal (e.g. HJD 2456127 and 2456149), 
and sometimes the sign of the signature is opposite (e.g. HDJ 
2455744, the 2456128-2456133 and the 2456176-2456177 se- 
ries). 

In order to search for possible significant variation of the 
LSD V profiles in the first data set, we have computed the stan- 
dard deviation of the individual V profiles. The results are plot- 
ted in Fig.|2](in black). If small but significant differences exist 
between the profiles of the sample, we would expect an increase 
of the standard deviation inside the profile with respect to the 
continuum where the signal is null. In Fig. [2] we observe that 
the standard deviation in the black profile is totally flat, which 
means that no significant variation is detected in the magnetic 
signature of HD 190073 over more than four years (2005-2009). 

The average and standard deviation of the V profiles of the 
second data set are overplotted in Fig. [2] (in red). In contrast to 
the first set, we observe an increase of the standard deviation 
inside the profile, well above the noise level, indicating that the 
profiles are not all similar and that the variations were definitely 
significant in 201 1-2012. 

In order to identify a possible modulation period of the ob- 
served variations in the 2011-2012 Stokes V profiles, we have 
computed periodograms by performing a least-squares fit of the 
B( curve with a sinusoidal function. A sine wave is expected 
if the magnetic field is similar to an inclined dipole. The peri- 
odogram computed with the 201 1-2012 B( values shows a clear 
minimum at 39.8 + 0.5 d with a reduced^ 2 of 3.6 (Fig. [3]). The 
resulting best sinusoidal curve is superimposed with the mea- 
sured B( values in Fig. [4] The best fit reproduces well all but one 
of the 2012 observations, confirming that a real modulation is 
present in our data set. However, a number of 2011 data point 
are not concordant with this period. We have therefore searched 
for a period in each of the 2011 and 2012 data sets separately. 
While in the 201 1 data no significant period could be identified, 
a period of 40 + 8 d is clearly detected in the 2012 set, with a re- 
duced^ 2 of 1.8, consistent with the value found using the whole 
201 1-2012 data set. 

The relatively satisfactory fit to a sine wave in 2012, while 
the 201 1 data fit such a wave quite poorly, together with the fact 
that the 2012 Bf data seem to vary strongly, with 8 out of 20 of 
the B( values below 0, while only 2 in 31 of the 2011 data are 
below zero, suggest that the variations observed during the two 
epochs are quite different. 




p (d) 



Fig. 3. Periodogram of the Be values measured in the Stokes V 
201 1-2012 profiles. The red horizontal line marks the 3<x signif- 
icance level above the minimum. 
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Fig. 4. Longitudinal magnetic field measurements (red squares) 
in the V profiles for the 2011 (top) and 2012 (bottom) data su- 
perimposed with the best fit (P = 39.8 ± 0.5 d, dashed lines). 



4. Discussion 

From 2005 to 2009 no variation of the Zeeman signature was de- 
tected. Catala et al. (2007) proposed three hypothesis: (i) the star 
is seen pole-on, (ii) the rotation and magnetic axes are aligned, or 
(iii) the rotation period is very long (longer than the time range 
of our observations, i.e. several years). No fundamental changes 
of the photospheric line profiles (in /) or in the circumstellar 
environment have been observed from 2005 to 2012. It is there- 
fore reasonable to assume that the rotation axis has not changed 
either. If the variations of the Zeeman signature are due to rota- 
tional modulation, as suggested by the periodic variations of the 
2012 data, hypothesis (i) can be ruled out. Moreover, we pro- 
pose the rotation period to be about 40 days and consequently 
hypothesis (iii) can also be ruled out. Therefore before 201 1, the 
magnetic axis was very likely aligned with the rotation axis (hy- 
pothesis (ii)). 

It seems reasonable to assume that the modulation of the 
Zeeman signature observed in 2011-2012 is at least partly 
caused by the rotation of the star. If the magnetic configuration 
at the surface of the star is stable over many rotation periods, a 
periodicity should be observed in the Zeeman signature. In 2012 
alone, all but one of the 20 Be data are well reproduced with 
a sine wave curve, implying that a periodicity is detected. We 
therefore assume that the field of HD 190073 in 2012 is roughly 
dipolar, and inclined with respect to the rotation axis, so that the 
variations are due mainly to oblique rotator-like rotation. In the 
2011 data, we are not able to find any periodicity using a sine 
wave, and the period that fits almost all the 2012 data does not 
fit the 2011 data at all well. The 2011 data, while not strongly 
variable, seem more chaotic than those from 2012. We suggest 
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Fig. 5. Magnetic (red squares) and non-magnetic (black points) 
Herbig Ae/Be stars plotted in a Hertzsprung-Russell diagram. 
The red square with error bars corresponds to HD 190073. 
The CESAM pre-MS evolutionary tracks for 1.5, 2, 3, 4, and 
6 M (black full lines), 0.01, 0.1, 1 and 10 Myr isochrones 
(blue thin dashed lines), and the zero-age main-sequence (black 
dot-dashed line) are also plotted. The birthline taken from 
Behrend & M aederl d2001l) is plotted with a blue thick dashed 
line. The limit over which the convective envelope is smaller 
than 1 % of the stellar mass, and the CC-birthline (see text) are 
overplotted with orange dot-dot-dot-dashed lines. 



that during the 2011 observing season the field structure was 
changing its intrinsic geometry, and variations due to rotation 
were only a minor part of the observed variation. 

The spectral lines of HD 190073 are very narrow, and their 
broadening may be dominated by turbulence, and only slightly 
affected by rotation (ICatala et al.1120071) . As a consequence we 
may have only a limited access to the rotational Doppler reso- 
lution inside the spectral lines. This Doppler resolution is useful 
to study the polarisation level per rotational velocity bin, which 
provides strong constraints on the magnetic configuration at dif- 
ferent longitudes and l atitudes (see details of Doppler imaging 
in e.g. lVogt et al.ll 1 987D . It is therefore an efficient way to detect 
complex fields. In the absence of rotational Doppler resolution 
we are unable to use this method, and both simple and complex 
fields can produce Zeeman signatures similar to those observed 
in HD 190073. 

We therefore propose the following scenario. Before 2009, 
the magnetic configuration of HD 190073 was stable and the 
magnetic axis was aligned with the rotation axis, resulting in 
an absence of modulation in the Zeeman signature. Between 
2009 and 201 1 a phenomenon occurred that perturbed the fossil 
field. As a consequence, the magnetic configuration lost its ax- 
isymetry, and a rotational modulation is now observed. In 201 1 
the magnetic field had not yet reached a stable dipolar configu- 
ration, while in 2012 the stability seems to have been achieved. 

As the photosphere and circumstellar environment of 
HD 190073 display similar characteristics throughout this pe- 
riod, there is no evidence that an external event, such as interac- 
tion with a companion having an eccentric orbit, or a collision 
with a massi ve comet (as observ ed sometimes in young Herbig 
stars, see e.g. lGrinin et al . 1996), occurred during that time, and 
perturbed the stellar surface field. The phenomenon that per- 
turbed the magnetic field of the star is therefore more likely of 
internal origin. 

We have investigated the predicted internal structure of the 
star, by comparing its position in the HR diagram with pre- 



MS C ESAM evolutionary track s (iMorell [19971) of solar abun- 
dance (lAcke & WaelkensTl2004l) . The distance of the system is 
not w ell known, because the Hipparcos parallax (Ivan Leeuwenl 
120071 0.70 + 0.74 mas) is very poo rly constrained. Th erefore 
we used the effective temperature of Catal a et al.l (120071 T e g = 
9250 + 250 K) and \ogg = 3.75 +0.3 to constrain the star's posi- 
tion in the HR diagram. The adopted value of log g was selected 
because our spectra can be reproduced equally well with log g 
of 3.5 or 4.0 when the effective temperature is varied within the 
error bars. In Fig. if we take into account the error bars, we 
observe that HD 190073 is situated just before or perhaps on 
the convective core birthline (CC-birthline), i.e. the condition in 
which the convective core has reached a mass larger than 1 % of 
the mass of the star. Before reaching the CC-birthline, the star is 
entirely radiative and hosts a fossil field. When the star reaches 
the CC-birthline, and the convective core starts to build u p, and 
convection is able to generate a dynamo magnetic field ( Mossl 



convection is able to generate a dynamo magnetic field (Moss 
1989; ICharbonneau & M acGregorl l200 U iBrun et al.ll2005l) . We 



expect this dynamo field to coup le with the fossil field at th e bor- 
der with the radiative envelope. iFeatherstone et all d2009) stud- 
ied, with 3D magneto-hydrodynamical simulations, the interac- 
tion between such a convective dynamo field and a stable ax- 
isymetric twisted dipolar fossil field in a surrounding radiative 
envelope. They showed that the coupling induces modifications 
of the fossil field geometry, in particular that the field becomes 
oblique. 

Therefore we propose that the dynamo field generated by the 
growing convective core has perturbed the fossil field anchored 
in the radiative zone of the star and is at the origin of the ob- 
served change in HD 190073. The simulations by Featherstone 
et al. (2009) included an already established dynamo in the core, 
perturbed by a fossil field in the envelope, while in HD 190073 
the fossil field is present first and the convective dynamo appears 
progressively. Therefore it would be useful to simulate this spe- 
cific case. In addition, the observed change occurred in less than 
2 years. This timescale is compatible with the proposed scenario, 
since stellar convective instabilities, and con sequently a dynamo, 
grow on timescales of about 1 month (e.g. lParker]|1975l) . 

To test whether the new magnetic configuration has already 
stabilised and to determine the precise parameters of this con- 
figuration, additional observations will be required. HD 190073 
represents a unique opportunity to witness the appearance of a 
convective core in a star. 
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